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Pionful nuclear effective field theory (NEFT) in the two-nucleon sector is examined from the 
Wilsonian renormalization group point of view. The pion exchange is cut off at the floating cutoff 
scale, A, with the short-distance part being represented as contact interactions in accordance with the 
general principle of renormalization. We derive the nonperturbative renormalization group equations 
in the leading order of the nonrelativistic approximation in the operator space including up to 0{p^) 
and find the nontrivial fixed points in the and ^Si-'^Di channels which are identified with those 
in the pionless NEFT. The scaling dimensions, which determine the power counting, of the contact 
interactions at the nontrivial fixed points are also identified with those in the pionless NEFT. We 
emphasize the importance of the separation of the pion exchange into the short-distance and the 
long-distance parts, since a part of the former is nonperturbative while the latter is perturbative. 



I. INTRODUCTION 

Nuclear effective field theory (NEFT) is the low-energy 
effective field theory [l| of nucleons based on symmetries 
of QCD, and is expected to give a modcl-indepcndcnt 
description of nuclear phenomena at low energies. Since 
the seminal papers of Weinberg 0-111, ^ ^ot of applications 
have been done with success for more than two decades. 
See Refs for reviews. It still has a wide range of 

phenomena to be explored. 

In spite of these successes, however, there are unset- 
tled issues at the fundamental level even in the simplest 
two-nuclcon sector; whether the pions should be treated 
perturbatively or not and how to treat the contact inter- 
actions. Because the understanding of the two-nucleon 
scattering is essential also for other phenomena, this issue 
is of central importance in NEFT. 

The original power counting due to Weinberg 
for the "effective potential" is nothing but the naive di- 
mensional analysis. The "effective potential" is plugged 
into the Lippmann-Schwinger equation and the scatter- 
ing amplitude is calculated, so that the pion exchange is 
iterated infinite times. In this Weinberg scheme the pions 
are thus treated nonpcrturbatively. Kaplan, Savage, and 
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Wise [l3l pointed out that the Weinberg power counting 
scheme is inconsistent, in the sense that the higher order 
counterterms are required to cancel the divergences at 
the given order. They (and, independently, van Kolck) 
proposed an alternative power counting scheme in which 
only the non-derivative contact interaction, the Cg-term, 
is treated nonpcrturbatively so that the pions are treated 
perturbatively [lTI - [l3l |. This scheme (widely known as 
KSW scheme) is free from the inconsistency, but Fleming, 
Mehen, and Stewart [l3| showed that the KSW scheme 
does not lead to the converging results in the channels in 
which the singular tensor part of the pion exchange con- 
tributes. Beane, Bedaque, Savage, and van Kolck [l^ 
proposed a remedy, known as the hybrid approach, in 
which the amplitudes are expanded around the chiral 
limit so that the 1/r^ part (which survives in the limit) 
is treated nonperturbatively. 

The treatment of the tensor part of the one-pion ex- 
change (OPE) is the source of controversy. See also 
Refs. |16l - [2l| for more details. 

Recently, Beane, Kaplan, and Vuorincn (BKV) [l^l 
considered a Pauli-Villars type regularization for the 
OPE in order to separate the short-distance part of 
the tensor interaction from the long-distance part. The 
short-distance part is represented by contact interactions. 
They employed the power divergence subtraction (PDS) 
renormalization and obtained the convergent results. 

Note that not all the \/r^ part of the OPE is singular. 
It is the short- distance part that is singular and spoils 
the convergence of the amplitude as found in Ref [l4| . It 
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means that all the l/r^ part of the OPE docs not need 
to be iterated. 

The idea of separating the singular short-distance part 
of the tensor interaction from the long-distance part and 
representing the former by contact interactions (BKV 
prescription) is essentially the Wilsonian renormalization 
group (RG) idea of effective interactions [l^, H^l- The 
separation scale in Ref. [2^ may be viewed as an analog 
of the floating cutoff in the Wilsonian RG. 

The Wilsonian RG is a useful tool to investigate the 
effects of the short-distance physics on the long-distance 
physics. When the quantum fluctuations with momenta 
higher than the floating cutoff A are integrated out, their 
effects are simulated by a series of local operators, which 
serve as low-momentum effective interactions. The cou- 
pling constants thus depend on the cutoff A, while the 
physical quantities such as scattering amplitudes do not. 

The Wilsonian RG has been applied to NEFT in order 
to examine the power counting issue [l^ [25l - f30j . (See 
also Ref. [3ll-l33l| for other use of the Wilsonian RG in 
NEFT.) The existence of the nontrivial fixed points of 
the RG equations (RGEs) in the S-waves explains un- 
naturally large scattering lengths. The power counting 
can be determined by the scaling dimensions of opera- 
tors around the fixed points. In particular, it is relevant 
operators that should be resummed to all orders. 

Birse [l^ examines the pionful NEFT by using the 
so-called "distorted- wave RG [2^" and claims that the 
scaling dimensions are shifted from those of the pionless 
theory due to the singularity of the tensor force in the 
spin-triplet channel. Also, his results imply that the ef- 
fects of pions do not decouple even at very low momenta. 

It sounds strange however that pions do not decouple 
at the momenta where the pionless NEFT is valid. The 
effects should be eventually represented by contact inter- 
actions at very low momenta. One should expect that the 
transition from the pionful NEFT to the pionless NEFT 
is smooth. A formulation of the Wilsonian RG which 
permits the smooth transition is desired. 

In this paper, we perform the Wilsonian RG analysis of 
the pionful NEFT in the two-nucleon sector, and exam- 
ine the effects of pions. In the Wilsonian RG approach, 
there is a single scale A which separates the long-distance 
physics from the short-distance physics so that one does 
not need to implement an additional regularization for 
the short-distance part of the OPE, as does in the BKV 
prescription. 

The key idea here is the separation of the short- 
distance part from the long-distance part. NEFT, as a 
low-energy effective theory, has a finite range of applica- 
bility, specified by the physical cutoff Aq. It is a general 
principle of EFT that the physics with the momentum 
scale larger than Aq is represented as local operators. 
The short-distance part of the OPE (S-OPE) should also 
be represented as local operators (contact interactions). 
(In the Wilsonian RG analysis, the floating cutoff A plays 
the role of Aq, after integrating the modes with momenta 
between A and Aq.) 



Note that, although the contact interactions are at- 
tributed to the effects of the heavy particles in the usual 
EFT lore, this is not precise. They arise also from the 
short-distance physics of light particles, such as pion 
exchanges. Note also that the representation of short- 
distance physics by local interactions is a general treat- 
ment and has nothing to do with how singular the S-OPE 
is. 

We find nontrivial fixed points in the and ^^i- 
■^Di channels, which are responsible for the large scatter- 
ing lengths. Importantly these fixed points are identified 
with those found in the Wilsonian RG analysis in the pi- 
onless NEFT. Thus the scaling dimensions at the fixed 
points are also the same. There is one relevant operator 
in each channel. That is, the pion interactions do not 
alter the scaling dimensions. 

We emphasize that the question of whether the pions 
are perturbative or nonperturbative is not well posed. 
Since the S-OPE is represented as contact interactions 
in the Wilsonian analysis, only the long-distance part 
of the OPE (L-OPE) is the proper interactions due to 
OPE. On the other hand, the S-OPE cannot be distin- 
guished from other contributions, such as heavier meson 
{p, uj, etc.) exchanges. (Because of the cutoff, we do 
not have enough resolution.) The RGEs tell us that, at 
the nontrivial fixed point, the L-OPE should be treated 
perturbatively, while there is a relevant operator (a part 
of which is the S-OPE) that should be resummed to all 
orders, i.e., nonperturbatively. 

The existence of the relevant operator at the nontrivial 
fixed point sharpens the distinction between the S-OPE 
and the L-OPE. If there were no relevant operator, the 
separation would not make much difference. 

Our Wilsonian RG permits the smooth transition from 
the pionful NEFT to the pionless NEFT. The nontrivial 
fixed points do not change nor the scaling dimensions. 
The L-OPE transmutes into contact interactions which 
represent the S-OPE as the cutoff is lowered. At the 
value of A lower than the pion mass, the most of the 
L-OPE has been changed into contact interactions, thus 
the system is well described by the pionless NEFT. 

The structure of the paper is the following. In Sec HI] 
we recapitulate the main points of the previous papers in 
order to introduce the notations and the main concepts 
in our analysis. In Sec. IIIII we discuss the nonrelativistic 
approximation. Starting with the nonrelativistic nucle- 
ons and the relativistic pions, we estimate the order of 
magnitude of a various kind of diagrams contributing to 
the RGEs and determine the leading terms. In Sec. IIVI 
wc present the RGEs and examine the structure of the 
fiows in the ^5*0 and "^Si-^Di channels. The nontrivial 
fixed points are found to be identified with those of the 
pionless theory. Sec |V] is devoted to the summary and 
the comments on the related works. In Appendix 1X1 the 
RGEs for the case of A < are presented and compared 
with the pionless case. In Appendix iBlwe discuss the sim- 
ilarity and the difference between the pionful NEFT and 
QED. 
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II. WILSONIAN RG ANALYSIS OF PIONLESS 
NEFT IN TWO-NUCLEON SECTOR 

In the previous papers [28l - [30| . we have explained the 
basic concepts of the Wilsonian RG and its relevance to 
the power counting in the NEFT. Here, we briefly recapit- 
ulate it, give some remarks, and introduce the notations 
used in later sections. 



A. What is the use of the Wilsonian RG in NEFT? 



nucleon system, the spin-singlet channel is considered 
to be in the weak-coupling phase, while the spin-triplet 
channel is considered to be in the strong-coupling phase, 
because of the (non)existence of a bound state in these 
channels. 

To summarize: the two-nucleon system with large S- 
wave scattering lengths is governed by the existence of 
nontrivial fixed points, and the Wilsonian RG is a sys- 
tematic tool to study the scaling dimensions on which 
the power counting should be based, around the nontriv- 
ial fixed points. 



The most basic idea behind the power counting is the 
order of magnitude estimate based on dimensional anal- 
ysis. In an EFT with the physical cutoff Aq, the dimen- 
sional analysis is usually based on this scale. At the clas- 
sical level, the (canonical) mass dimension of an operator 
is determined with respect to the kinetic term. For exam- 
ple, a Dirac field ip with a kinetic term Ckin = i)i$'4> has 
mass dimension three halves, [ip] = 3/2. The dimensions 
of other operators are determined accordingly. The op- 
erator (ipip)'^ has dimension six, so that it enters in the 
Lagrangian as {c/ Aq){iPiP)'^ , where c is a dimensionless 
constant. The coupling constant c/Aq associated with 
the operator (jpip)'^ has dimension —2, which counts the 
power of Aq. 

Quantum fluctuations may, in general, change the clas- 
sical dimensional analysis. The quantum counterpart of 
the (canonical) dimension is called the scaling dimension, 
which can be obtained by the RG analysis. Wilsonian RG 
is a nonperturbative tool to handle the quantum fluctu- 
ations. 

An operator whose coupling has a negative (scaling) 
dimension is called irrelevant because it becomes less 
important at lower energies. An operator whose cou- 
pling has a positive (scaling) dimension is called relevant 
because it becomes more important at lower energies. 
An operator whose coupling is dimensionless is called 
marginal. The scaling dimension of an operator is the 
measure of how important it is. It is therefore natural to 
consider the power counting on the basis of the scaling 
dimensions. 

In the S-wave scattering of two nucleons, the scattering 
lengths are known to be much larger than the "natural" 
size, I/Aq. (In the case of the pionless NEFT, the phys- 
ical cutoff is of order of the pion mass, Aq ^ ©(vtItt).) 
From the RG point of view, the "fine-tuning" is related 
to the existence of a nontrivial fixed point (and a critical 
surface) of the RG flow. 

Around the nontrivial fixed points the scaling dimen- 
sions are drastically different from the canonical dimen- 
sions. It has been shown [2^, |2^, Is^l that the coupling 
which corresponds to the scattering length becomes rel- 
evant, although it is irrelevant at the classical level. 

There are values of coupling constants with which the 
scattering length is infinite. This set of coupling con- 
stants forms a critical surface. It separates the weak- 
coupling and the strong-coupling phases. For the two- 



B. Remarks on Wilsonian RGBs with Galilean 



invariance 



There are several formulations for the Wilsonian RG 
analysis [s^ - fsoj . which are however esscntailly equiva- 
lent. The most popular one is the functional RG method. 
(See Refs. [iO,[iH for reviews.) In this formulation, a cut- 
off function is introduced for each propagator to suppress 
the low-frequency fluctuations. The effective averaged 
action Fa , which interpolates the bare classical action 
S'[$] (A = Ao) and the effective action r[$] (A = 0), 
depends on the floating cutoff scale, A, as 



dTA 
dA 



1. 



-Tr 



dRA 

-dA ^^(2) 



RaV' 



(2.1) 



where $ is the classical field, F(2) stands for the sec- 
ond derivative of the averaged action Fa with respect to 
<f>, and R\ is the cutoff function which suppresses the 
fluctuations with p < A. Note that, although it is an 
"one-loop" equation, it contains all the nonperturbative 
information. 

A straightforward application of this formulation to a 
nonrclativistic system however encounters difficulties. In 
the usual formulation for a rclativistic system, one con- 
siders the theory in Euclidean space. The cutoff is im- 
posed on the magnitude of four-momentum of the prop- 
agator. On the other hand, in a nonrclativistic system, 
space and time should be treated differently, and thus the 
Euclidean formulation cannot be used. One might rather 
like to consider a cutoff imposed on the three-momentum 
in the propagator. But such a cutoff necessarily breaks 
Galilean invariance of the nonrclativistic system. There 
is no obvious way to impose a Galilean invariant cutoff 
at the averaged action level. This is a general feature 
independent of the choice of the cutoff function. Fur- 
thermore, if the cutoff function is not smooth enough, 
non-analytic terms in momenta arise0. See Ref. [l^ for 
an example in a similar context. 



^ If one considers only a few leading order terms in the derivative 
expansion, typically in the local potential approximation, non- 
smooth cutoff does not cause the problem. It is the reason why 
the problem of non-analyticity is not revealed in many applica- 
tions. 
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In a system of two nonrelativistic particles, there is a 
simple and physically sensible way out of this problem: 
it is to impose a cutoff on the relative three-momentum 
of the two particles, which is Galilean invariant. In ad- 
dition, the results are very insensitive to the choice of 
the cutoff function. See Appendices of Refs 0, [s^l ■ In 
particular, the results with a sharp cutoff are the same 
as those with a smooth one. It is a technical advantage 
that a sharp cutoff can be used because it simplifies the 
calculations considerably. 



Pionless NEFT up to ©(p^) in the ^So and 
^Si-^Di channels 



In Rcf. [29|, we consider the pionless NEFT without 
isospin breaking. The relevant degrees of freedom are 

I 



nonrelativistic nucleons, which interact with themselves 
only through contact interactions. In the two-nucleon 
sector, they are four-nucleon operators with an arbitrary 
number of derivatives. An operator with more derivatives 
has higher canonical dimensions than the ones with less 
derivatives. 



Since there arc infinitely many operators involved in 
the fiow equation, one needs to introduce a truncation of 
the space of operators in the averaged action in order to 
solve it. We retain only the operators with derivatives up 
to a certain order. We simply count the number of spa- 
tial derivatives (V ^ p) and a time derivative is counted 
as two spatial derivatives {dt ^ p^). We consider the 
following ansatz for the averaged action up to 0{p'^), 



(T).o(T) 



2S(S)0(SB) 



(^S'o channel) 
{^Si-^Di channel) 



where the operators in the ^S'o are given by 



of) 



[SB) 



N^P^,^^ {idt + —\N\ ( A^Pi^'A) + h.c 



and in the '^S'l-'^Di channel. 



O, 



(T) 



(T) 



{n'^p^^^n)\n'^p'>^^n) 



n'^pPn] [N' P, 



(SD) 



(TB) 



)\N^PP'^^N)+h.c 

/^^^)a)^{a-(v.v,4..,v^)p 

^ A^l (a^'^P^^^^a) + h.c. 



n^pP 



N} +h.c. 



idt + 



(2.2) 



(2.3a) 
(2.3b) 

(2.3c) 

(2.4a) 
(2.4b) 

(2.4c) 

(2.4d) 
(2.4e) 



wfierewe have introduced the notation V ^ = + V — 
2v • V and the projection operators [Tij. 

piS) ^ 1 pP ^ ^a^a\\ (2.5) 

for the ^S'o (spin singlet) channel and the ^Si (spin 



triplet) channel respectively. The nucleon field N{x) 
with mass M is an isospin doublet nonrelativistic two- 
component spinor. Pauli matrices cr' and act on spin 
indices and isospin indices respectively. The two chan- 
nels are completely decoupled, and thus we can consider 
each channel separately. 
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Note that the possible forms of the operators are re- 
stricted by Galilean invariance. Note also that we have 
included "redundant operators," of^^ and o'^^\ be- 
cause they are necessary in a consistent Wilsonian RG 
analysis j28| . 

We introduce the following dimensionless coupling con- 
stants, 



MA is) 



MA3 



(S) 



A3_ 
27r2 



B^^\ (2.6) 



for the spin-singlet channel and 



MA 



2^^o 
MA3 4 



(T) 



MA'^ 



(T) 



A3 

27r2 ' 



27r2 3 



:c: 



iSD) 



(2.7) 



for the spin-triplet channel. With the sharp cutoff on the 
relative momenta, the flow equations that determine the 
dependence on t = ln(Ao/A) of the coupling constants 
can be written as 1301 



—j--+dcxc = y xAXBi^ ~— — — 



c 
(2-8) 

where xc stands for one of the dimensionless coupling 
constants, and dc is the power of A in the definition of 
the dimensionless coupling constant. In the following, we 
call —dc the canonical dimension0of the coupling. A{P) 
is defined as 

^^^^ = " S"' ^(^^ = MA{P)/A^ (2.9) 

where P = (P", P) is the total momentum of the system. 
FA{Pi,Pf) is the momentum-dependent factor associated 
with the coupling xa- 



F. 



F, 



-27r2 
MA ' 



A3 Z^"^*' 



-27r2 



4AfA3 



(ri2 + r34) ■ 



(2.10) 



1=1 



for the spin-singlet channel, and 
-2^2 



A/A 



-2^2 



^■34) , 



F,- 



^3 Z^*^'' 

i=l 

3 /-2^2x 

4 1^ MA3 j 



K2 P{2 + P34 P34 - o^'^ (^12 + ''34) 



(2.11) 



This definition of the canonical dimension reflects the nonrela- 
tivistic scaling property. See Ref . [42l . 



for the spin-triplet channel, with 



(p.,-PjY. (2.12) 



Pi and P2 are the incoming momenta of the nucleons to 
the vertex, and pa and p4 are the outgoing momenta from 
the vertex. The notation \c stands for the operation of 
taking the coefficient of Fc(j>i,Pf) in the sum. 

By using the formula Eq. (|2.8p the RGEs are obtained. 
The explicit expressions are not presented here, but can 
be read easily from the RGEs (Eqs. (|4.6p and Eqs. (|4.12p ) 
for the pionful theory discussed in Sec. IIVI 

There is a nontrivial fixed point in each channel, which 
is relevant to the physical two-nuclcon system: 



ix\y\z'') = 



1 1 

""^'^2' 2 



in the spin-singlet channel, and 



(2.13) 



(2.14) 



in the spin-triplet channel. At the nontrivial fixed point, 
the operators get large anomalous dimensions and there 
is one operator that becomes relevant with the scaling di- 
mension of the coupling constant being one in each chan- 
nel. 



III. NONRELATIVISTIC APPROXIMATION, IR 
ENHANCEMENT AND THE LEADING ORDER 

IN A/M 

In this section, we consider the inclusion of pions as 
dynamical degrees of freedom. It extends the range of 
applicability of NEFT to higher momenta beyond the 
pion mass scale. The contact interactions in the pionless 
theory resolve into the effects by pion propagation and 
the rest. The physical cutoff Aq is now larger than m-^, 
and we suppose that Aq is of order 400 MeV. 



A. Chiral symmetry and nonrelativistic nucleons 

The most important feature of the pionful theory is 
chiral symmetry, SU{2)l x SUiTjn spontaneously broken 
to SU{2)v- It is convenient to introduce the field E, 
which transforms linearly as 



(3.1) 



where L and R are the elements of SU{2)l and SU{2)ii 
respectively. The pion field 7r"(x) may be defined through 



S(x) =exp(*7r'^(a;)T7/), 



(3.2) 



where / is the pion decay constant in the chiral limit. 
The nucleon field transforms as 



N{x) U{x)N{x), 



(3.3) 
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where U{x) is a function of L, R, and E(x) and defined 
through 

ax) L^xWixy = Uix)ax)R\ (3.4) 

where £,^{x) = T.{x), i.e., £,{x) = exp(i7r''(x)T'^/2/). 

The chiral invariant Lagrangian for the nonrelativistic 
nucleon interacting with pions is given as 



-NR 



iDo 



2M 



N + gAN'fcT ■ AN 



-I?W!^Tr(I] + St)oW 



where the chiral covariant derivative is defined as 



(3.5) 

IS 

(3.6) 



and and are defined as 



V, 



M — r, 

i 



(3.7) 
(3.8) 



The superscript (c) denotes the channel to be speci- 
fied, and the summation is also implied, i.e., the term 
02'^^ 0^2^ for the spin-triplet channel contains cl^^^O^'^ 

and C^^^0^2^\ Note that the derivatives in the four- 
nucleon operators should be replaced by the covariant 
derivatives, though we do not explicitly show them here. 

The ellipsis in Eq. (|3.5p denotes other terms, e.g., the 
higher order operators, including six-nucleon operators, 
four-nucleon operators with more than two derivatives, 
etc. It also contains the counterterms of the form N'^ N 
and N^W'^N, which play an important role in the renor- 
malization of the nucleon mass. Chiral symmetry does 
not prevent the appearance of the operators that con- 
tain more than one A, such as N'^A^N. But it is not 
generated by nonrelativistic reduction from the simple 
(chiral invariant) Dirac action. It implies that the oper- 
ators have smaller coefficients than 1/M. Furthermore, 
we will see that it either does not contribute to the RGEs 
for the four-nucleon operators, or gives only suppressed 
contributions to the order we are working. 

The Lagrangian for the pions is given by 

= ^ [Tr (a^S^a'^E) + m^Tr (S^ + e)] + • • • . (3.9) 



B. Order of magnitude estimation of the 
contributions to tlie Wilsonian RGEs 

In order to perform the Wilsonian RG analysis, one 
usually needs to use a cutoff function that preserves all 
the symmetries of the theory. Unfortunately, there does 



not seem to exist a manifestly chiral invariant cutoff func- 
tion which controls all the fluctuations, because of the 
nonlinearity of the transformation (43l |. Furthermore, it 
is known that perturbation theory generates apparently 
noninvariant terms (ANTs) even with the lattice and the 
dimensional regularization, which preserve chiral sym- 
metry [13, ■ ANTs are also expected to appear in the 
Wilsonian RG analysis, but it is not obvious how to treat 
them. 

In addition, the inclusion of pions makes the notion of 
relative momentum obscure. Furthermore, because pions 
are relativistic, Galilean invariance does not make good 
sense as a constraint. 

Fortunately, however, it turns out that these problems 
do not interfere with the leading order calculations in the 
nonrelativistic approximation. The argument is based on 
the order of magnitude estimate of the contribution of 
each diagram to the Wilsonian RGEs. 

We are going to obtain the RGEs for the two-nucleon 
sector in the next-to-leading order in momentum expan- 
sion (to 0{p^)) of the averaged action in the nonrelativis- 
tic formulation. Let us first remind the following things: 

• The contributions to the Wilsonian RGEs only 
come from one-loop diagrams. 

• Because of the nonrelativistic feature, diagrams 
with anti-nucleon lines are absent. Thus, the di- 
agrams arc divided into sectors, each of them is 
specified by the nucleon number. The rt-nucleon 
sector has n nucleons at each time slice. 

• The contributions to the two-nucleon sector do not 
come from 71-nucleon operators with n > 6. 

• Our primary concern is the rcnormalization of the 
four-nucleon operators with no pion emission and 
absorption. Chiral symmetry constrains the rcnor- 
malization of the operators in which the pion fields 
appear through covariant derivatives. 

• Supplementarily, one needs to consider the self- 
energy diagram of the nucleon and the diagrams 
for the nucleon-pion vertex that contribute to the 
RGEs for the four-nucleon operators. 

Note also that there is no contributions to the pion 
self-energy to this order. 



1. Four-nucleon operators 

According to the above-mentioned remarks, we need to 
consider the diagrams given in Fig.[T]for the four-nucleon 
operators. 

In order to obtain the RGEs, we need to evaluate the 
contributions from the so-called "shell-mode," the loop- 
contributions with the magnitude of the loop (relative) 
three- momentum k = \k\ is between A — dA and A. The 
argument is similar to the "power counting" with the 
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FIG. 1. Contributions to the four-nucleon operators. Tiie 
four-nucleon vertices represent contact interactions collec- 
tively. The dotted lines represent the pion propagators. For 
the diagrams (e), (f), (g), (h), and (1), there are also mirrored 
diagrams with the left and the right interchanged. 



scale Q found in the literature, but here, (i) we do not 
consider the amplitude but the contribution to the RGEs, 
and (ii) the magnitude of the momentum in the loop is 
actually A, while in the "Q-counting" it is assumed that 
the dominant contributions come from the loop momen- 
tum of order Q. 

To see how the dominant contributions to the RGEs 
emerge in certain diagrams, let us consider the diagram 
(b) in Fig. [T] as an illustrative example. At the moment, 
we assume that the floating cutoff A is larger than 111^^ 
for simplicity. The diagram contains the loop integral 



hell (27r)4 fc2 - m2 ze (pO + _ (p + fc)2/2Af + ; 



- fc") - (p' - fc)2 /2M + ie 



(3.10) 



where /g^cn denotes the shell mode integral with the 
restriction on the magnitude of the relative three- 
momentum. There are four poles in the complex fc" 
plane, one nucleon and one pion poles in both the upper 
and the lower half planes. One can evaluate the integral 
by the residues of the poles in either half plane, and finds 



that the nucleon pole gives a dominant contributions, 
d^k -fe2 



1 



{p°+p'°) - (p+fe)72M- (p'-fc)2/2M 
d^k -{k- {p^p')l2f 



hell (2^)3 [_pO + (fe + p/2Y/2Mf - ujI 



1 



E - P^/AM -k'^/M 



k~(p-p')/2 

(3.11) 



where Wfc = \/k + , and E = +p'" and P = p + p' 
are the total energy and the total three-momentum of 
the two nucleons, respectively. In going to the second 
line, we have made a shift of the integration momentum 
so that k is now the relative momentum. Since M 3> 
|fc| = A > \P\, p°, and E - P^/iM < A^/Af , one may 
estimate the loop integral as 

-1 



27r2 



-MdA. 



(3.12) 



The pion pole gives ^ AdA/Air'^, which is smaller than 
the dominant contribution by a factor of A/M. 

On the other hand, the diagram (e) in Fig. [T] contains 



d^k 



ik^ 



shell 



(27r)4 fc2 - m2 + ie {p° - /fc") - (p - fc)2/2A/ + ie 



{p'° - fcO) - (p' - fc)2/2Af + ie ' 



(3.13) 



which is similar to Eq. p.lOp . but there is a crucial 
difference. No nucleon pole appears in the lower half 
plane. Thus the contribution can be evaluated by the 
residue of the pion pole in the lower half plane and gives 
^ — AdA/47r2. This is suppressed by a factor of K/M 
compared with the diagram (b) in Fig. [TJ 

It is a general feature that the two-nucleon reducible 
(2NR) diagram acquires the so-called "IR enhancement," 
first noted by Weinberg 0] , where the residue of the nu- 
cleon pole gets an enhancement factor A// A. Note that 
the IR enhancement is a consequence of the nonrelativis- 
tic kinematics. In addition, with the residue of the nu- 
cleon pole, the pion propagator becomes 

(3.14) 



[-p° + {p + kY/2MY -ujI 



Noting p° - p2/2A/ < A2/A/ < A, it is approximated 

by 



(3.15) 



fe2 



for fc ^ A. That is, the effects of the pion can be repre- 
sented as the instantaneous potential. 

Although we have shown that the mechanism of "IR 
enhancement" works in the case of < A, it is actually 
independent of this assumption, and it also works in the 
case of A < m^. 



8 



(a) 



(b) 



(c) (d) 
FIG. 2. Contributions to the ^A-vertex. 



FIG. 3. Contributions to the self-energy of nucleon. 



2. Vertex corrections and the nucleon self-energy 

We also need to consider the renormalization of the 
g^i-term. The relevant diagrams are depicted in Fig. [5] 
There are contributions to the self-energy of the nucleon, 
shown in Fig. [S] These are potentially important to the 
renormalization of the four-nuclcon operators. 

The examples in the previous subsection imply that 
the nucleon self-energy diagram and the contributions to 
the nucleon-pion vertex do not have IR enhancement. We 
will see shortly that these contributions can be neglected 
in our leading order calculations. The key point is that 
the appropriate dimensionless coupling constant for the 
nucleon-pion coupling is, as wc will explain in the next 
section (Sec. IIV A[) . given by 7, 



_ MA / gA 
27r2 ^^2/ 



(3.16) 



The contributions from the nucleon-pion vertex to 7 
should be written in terms of 7. 

There are four diagrams contributing to the nucleon- 
pion vertex shown in Scc. lIIIB] The tadpole diagram (b) 
can be absorbed in the redefinition of the coupling con- 
stant gA- Each of the diagrams (c) and (d) has an extra 
1/M factor because the two-pion vertex comes from the 
kinetic term £)^/2M so that they are suppressed com- 
pared to the diagram (a). The diagram (a) gives rise to 
the second term of the right-hand side of the following 
RGE expressed in terms of 7, 



dt ' \M ' ' 



(3.17) 



(The contributions from the nucleon wavefunction renor- 
malization is not considered here.) Note that the last 
term has an explicit A-dependence. In Eq. (|4.7p the last 
term has been neglected because A/M is much smaller 
than one. 

The self-energy diagram S(p) shown in Fig[3]gives con- 
tributions 



A 



dA 



A + Bp" + Cp^ + OiA^/Ad^), (3.18) 



where A, B, and C are the constants that depend on the 
couplings and the cutoff A. The constants A and C are 
canceled by the counterterms, leaving the pole of the nu- 
cleon propagator intact. After doing so, the propagator 
with this shell-mode contribution becomes 

-1 

(3.19) 



i-b'4 

A 



P 



2M 



so that the contributions to the wave function renormal- 
ization constant for the nucleon field, Zn, can be written 
as 

dZN 



dt 



= B. 



The order of magnitude estimate of B gives 

A_ 



B 



7, 



(3.20) 



(3.21) 



so that the inclusion of the effect of wavefunction renor- 
malization does not alter the results of the leading order 
calculations. 



C. Averaged action with L-OPE 

We have seen that the leading contributions to the 
RGEs consist of the one-loop 2NR diagrams with contact 
interactions and/or the instantaneous pion exchanges. 
Actually these contributions can be generated by a much 
simpler action than that in Eq. p.Sp . Thus we start with 
the following ansatz for the averaged action Fa, 



d^x^^-Di^^m 



9\ 



^ I dt I d^'x d^y 



0^''\x,y)Vl 



4/ 

+0lr(x,y)V 
-60iJ\x,y) (d^d^^h^Wl 



where we have introduced 



Y{\x-y\), 

(3.22) 



O(^) {x, y) = (iV^(x)Pf) iV(y)) (y) P^'^'' N [x] 



(3.23) 



0^;f\x,y)^(N^{x)PPN{y))\N^{y)PpN{x) 



(3.24) 
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and 



Y(r) 



1 e- 



(3.25) 



Note that all the derivatives are now the usual ones, not 
the covariant derivatives. It means that chiral symme- 
try is broken explicitly, but the breaking is of higher or- 
der in A/M, as is seen from the derivation. Note also 

that the operator corresponds to Dj^-* is the same as 

(c) 

that to Cq , but the former is a part of the operator 
that emits/absorbs pions and is of higher order in the 
p-expansion than the latter. 

The effects of pions are represented as one-pion ex- 
change interactions. It is the L-OPE because the av- 
eraged action is defined with the cutoff A, so that the 
S-OPE (with the momenta larger than A) is included in 
the contact interactions. 

The last term in Eq. (j3.22[) represents the tensor force 
of L-OPE in the spin-triplet channel. 



IV. WILSONIAN RGES FOR PIONFUL NEFT 

0(p2) IN THE \So AND '■^Si-^Di CHANNELS 

We are now ready to calculate the RGEs for the two- 
nucleon system in the S-waves. The formula Eq. (|2.8p 
can be used with slight modification. 

• We introduce the dimensionless coupling constants 
7 defined in Eq. p.l6p for L-OPE, and u and u' for 



defined by 



27r2 



A/A3 



-D. 



(T) 



(4.1) 



Note that qa and / appear in Eq. p.22p only 
through the combination of [gA/ fY- 

• The momentum-dependent factors associated with 
7, u and u' are given by 



Fis{Pf,Pi) = 



/-27r2^ 




V MA ) 





'"13 



ri3-fm^ 



ri4-|-m^ 





\ 




V MA ) 




ri3 + ml 



ri4 + ml 



F — F 



'2iT'^ml 
MA^ 



(4.2) 

(4.3) 
(4.4) 



Since the F^s and F^p have a bit more compli- 
cated momentum dependence than those in the 
pionless theory, the formula contains nontrivial 
integrations over angular variables. The part 



FA{Pr,A)FB{A,pf) in Eq. 
{FA{pt,A)FB{A,pf)), where 



is replaced by 
is defined as 



1 

47r 



(4.5) 



where stands for the angular variables of fc. See 
the original derivation of the formula for the de- 
tails j3G |. 

In the following, we assume that m^^ < A and expand 
in powers of mT, / A. The case of A < 771,^ is discussed in 
Appendix [X] 



A. Spin-singlet channel 

In the spin-singlet channel, we have the following 
RGEs: 



dx 
It 



dy 

dt 



+ 2a;?; + + 2xz + 2yz + 



-2(x -I- ?/ -I- z)7 — 7^ 
-32/ 



^.T^ + 2xy + -^yz- 



-(x + 2?/)7- ^7^, 



(4.6a) 



(4.6b) 



dz 

— = -3z 
dt 



1 'J 1 9 So 

—x" + XV H V — xz — vz z 

2 ^2 2 



1 



+ (.T + y-z)7+-7 



(4.6c) 



^ = -3u -2{x + y + z)(u - 7) - 2u7 + 27^ (4.6d) 
dt 



and for 7, 



c?7 



-7- 



(4.7) 



The first lines of Eqs (|4.6ap - (|4.6cp arc the same as 
those in the pionless calculations obtained in Ref. [29j . 
The terms in the second lines express how the L-OPE 
is rearranged into the S-OPE when the floating cutoff is 
lowered. 

We emphasize the choice made here of the dimension- 
less coupling constant 7. There are several ways to make 
a dimensionless quantity from the combination {gA/ fY', 
M , and A. Our choice is the one for which the RGEs for 
the coupling constants of the contact interactions do not 
have explicit A-dependence. If the explicit A-dependence 
were present in the RGEs, the iterative property (self- 
similarity) would be lost and the concept of fixed points 
would become obscure. Since the unnaturally large scat- 
tering lengths in the S-wave scattering are believed to be 
related to the nontrivial fixed points, one needs to use 
such dimensionless variables that allow fixed points. 
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The nontrivial fixed point of tlie RGEs (|4.6p and 
Eq. (j4.7p relevant to the real two-nucleon system is found 
to be 



(x ,?/ ,z ,u ,7 ) 



1 1 



(4.8) 



which is identified with that found in the pionless NEFT, 
given in Eq. ^J^. 

We linearize the RGEs around the fixed point, and 
find the following eigenvalues (scaling dimensions) and 
the corresponding eigenvectors: 



1^1 



+1 : ui 



V3, 



Us = 



-1 



V / 

-2 


V J 



V2 



-1 : U2 



-1 : U4 = 



1 


V / 

* 



1 

Vo/ 



(4.9) 



The eigenvalues vi, V2, and and corresponding eigen- 
vectors Ml , U2, and Us can be identified with those found 
in the pionless theory. Therefore, the power counting is 
not modified by the inclusion of the pions; only the one 
relevant operator iui) should be resummed to all orders. 

Note that the eigenvalue problem is five dimensional 
and one thus expects five pairs of eigenvalues and eigen- 

I 



vectors. It is easily seen that the eigenvalue —1 is triply 
degenerate. Two of the (nonzero) eigenvectors are U2 and 
U4, but the third one does not exist. This is not a math- 
ematical inconsistency, however. Although we do not 
understand very well the reason why the third eigenvec- 
tor is missing, it clearly has something to do with the 7- 
direction, because the vector in the 7-direction cannot be 
expressed as a linear combination of u^'s (i = 1, • • • , 4). 

Eq. (|4.7p shows that the L-OPE is irrelevant. It im- 
plies that the L-OPE should be treated as a perturbation. 
Note that there is a typical scale in the pionful NEFT, 

Aatjv, 



J-^NN — — I 



M \9A 
and our 7 is related to it as 



7(A) = 



A 



A 



NN 



(4.10) 



(4.11) 



Kaplan, Savage, and Wise [TT!, [Ill regard p/Aattv as an 
expansion parameter. Our finding is consistent with their 
approach. 



B. Spin-triplet channel 

In the spin-triplet channel, we have the following 
RGEs: 



~dt 



dt 



dw' 
~dt 



du' 
~dt 



x'"^ + 2x'y' + y'^ + 2x' z' + 2y' z' + z'^ + 2w'^ 



-2{x' + y' + z' - 4u7')7 - 87^ 



1 ,2 



2xV + ^/ 



V'z' 



1/2 ,2 

-z +w 



-(x' + 2y')7+^7', 



^ = -3z' 



Ix^' + xV + ^y'' 



x' z' — y z' 



+ {x' + y' - z' - 4u;')7 + 27^ 



-3w' 



1 



X w + y w + z w 



+ 7(2x' + 2y' + 2z' - 9it;')7 + 
5 

-3ii' - 2{x' + y' + z')v! + 2{x' + y' + z' - 4u;')7 - 2u'j + 187^. 

I 



(4.12a) 
(4.12b) 

(4.12c) 

(4.12d) 
(4.12c) 



Here, again, the first lines of Eqs (|4.12ap - (|4.12d|) are in Ref. [2^, while the second lines are the contributions 
the same as those in the pionless calculations, obtained 
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from L-OPE. 

Note that the magnitude of the coefficients of 7^ is 
large compared to those for the spin-singlet channel. This 
is the effect of the tensor part of L-OPE. 

The nontrivial fixed point of the RGBs (j4.12[) and 
Eq. (j4.7p relevant to the real two-nucleon system is found 
to be 



-1, 



1^,0,0,0 



(4.13) 



which is the same as that found in the pionless NEFT, 
given in Eq. (|2.14|) . 

The RGEs linearized around the fixed point lead to the 
following set of eigenvalues and the eigenvectors: 



i/i — +1 : ui = 



Vs. 



U3 



V5 



1 : M5 



/ 1 \ 

1 

-1 



V J 

( 2 \ 

-1 

-2 



V ) 







V2 



-1 : U2 = 



1/4 



-2 : U4 



/ \ 

-1 

1 



V / 





1 



Vo/ 



(4.14) 



Here, again, the eigenvalues 1^2, 1^3 and and corre- 
sponding eigenvectors wi, U2, U3 and M4 can be identified 
with those found in the pionless theory. Therefore, in 
spite of the presence of the tensor force in this channel, 
the power counting is not modified. As in the spin-singlet 
channel, only the one relevant (ui) should be resummed 
to all orders. 

As in the spin-singlet channel, the eigenvalue —1 is 
triply degenerate, but one eigenvector is missing. 



C. The effects of pions 

In this section, we discuss the effects of pions to the 
RGEs and thus to the power counting. 

As we discussed in Sec. Ill G[ the scaling dimensions at 
the nontrivial fixed point determine the power counting 
of NEFT in the S-wavcs in the two-nucleon sector. In 
the previous sections, we have seen that the nontrivial 
fixed points as well as the scaling dimensions at them in 
both ^5*0 and '^S'l-^Di channels remain the same even 
when pions are included as explicit degrees of freedom. 
Therefore the power counting for the contact operators 



in the pionful theory is the same as that of the pionless 
theory. 

The L-OPE does not affect the location of the fixed 
points and the scaling dimensions there, but the details 
of the flows. The "strong" tensor force in the spin-triplet 
channel does not change the properties of the nontriv- 
ial fixed point, but gives rise to the strong 7-dependence 
(because of the large factor "6" in Eq. (|3.22p ) in the chan- 
nel and modifies the flow considerably compared to the 
spin-singlet channel. 

The effects can be read off from the 7-dependence of 
the RGEs. To see the effects in more details, let us in- 
troduce new variables. 



X = X + y + Y = y + Z = y 



(4.15) 



and rewrite the RGEs in Eq. (j4.6p for the spin-singlet 
channel as follows: 



dX 
~dt 
dY 
dt 
dZ 
lit 

du 



-X-2Y-XY-X'^-{2X + Y)-f--/^, (4.16a) 

-{3 + X + j)Y, (4.16b) 

-3Z - X{X -Y + 2Z) - {2X - Y + 2Z)j - 7^ 

(4.16c) 

-3u-2Xu-2{X + u)-/ + 2-f'^. (4.16d) 



We see that the RGEs for {X, Y, 7) form a closed subset 
and can be solved without knowing the ffow for Z and u. 
In these variables, the nontrivial fixed point is given by 

(x^r^7'^) = (-1,0,0). 

Similarly, we introduce 

X' = x' + y' + z, Y' = y' + z', Z' = y' - z' , {4,.n) 

in the spin-triplet channel, and the RGEs can be rewrit- 
ten as 



dX' 
IT 

dY' 

~dr 

dZ' 
~dt 

dw' 
~dt 

(hj!_ 
dt 



-X' - 2Y' - X'Y' - X'"^ - 2w'^ 

-(2X' + r'-4^i;')7-7^ (4.18a) 

-(3 + X')Y' - [Y' + 4u/)7 + 87^, (4.18b) 

-3Z' - X'{X' -Y' + 2Z') - 2w'^ 
-(2X'-y + 2Z'-4w')7-7^ (4.18c) 



~3i/;' — X'w' 



i(2X'-W)7 
5 



+272, (4.18d) 
-3u' - 2Xu' + 2{X' -u' - 4w')7 
+I872. (4.18e) 



The RGEs for (X',F',w',7) form a closed subset, and 
the nontrivial fixed point is given by {X'* , Y'* , w'* , 7*) = 
(-1,0,0,0). 
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FIG. 4. (Color online) The surfaces which separate the ini- 
tial values into those from which the flows go to the strong 
coupling phase (the left side), and those to the weak coupling 
phase (the right side) . The surface on the right is for the spin- 
triplet channel, and the other on the left for the spin-singlet 
channel. All the points have the initial values with w' = for 
the spin-triplet channel. 



Thanks to the similarity of the RGEs in both channels, 
we can project the flows onto the {X, Y, 7) and {X', 1^' , 7) 
space and compare them in the single plot. In Fig. 21 we 
have drawn the surfaces which separate the initial values 
into two regions: the region from which the flows go to 
the strong coupling phase, and the other to the weak 
coupling phase. It is evident that the region of initial 
values to the strong coupling phase is larger in the spin- 
triplet channel than in the spin-singlet channel. 




FIG. 5. (Color online) Typical RG flow lines in the spin- 
singlet channel. The intervals between points indicate how 
fast they flow. 



In Figs. [5] and [HI we show some typical flows. Note 
that the large coefficient of the 7^ term of Eq. (j4.18bp , in 
comparison with Eq. (|4.16bp results in the large bending 
of the flow lines in the Y' direction. 




FIG. 6. (Color online) Typical RG flow lines in the spin-triplet 
channel with the initial value of w' set equal to zero. The 
intervals between points indicate how fast they flow. Note 
that the strong "dragging force" in the Y' direction. 



V. SUMMARY AND DISCUSSIONS 
A. Summary 

In this paper, we consider the pionful NEFT in the 
two-nucleon sector in the S-waves in the leading order 
of the nonrelativistic approximation in order to study 
the power counting issue from the Wilsonian RG point 
of view. We show that the leading order contributions 
to the RGEs come from the two-nucleon-reducible dia- 
grams, and the pion propagators are dominated by the 
instantaneous Yukawa potential. 

The separation of the pion exchange contributions into 
the L-OPE and the S-OPE is emphasized on the basis of 
the general effective held theory philosophy. The L-OPE 
is expressed as the Yukawa potential in the averaged ac- 
tion, while the S-OPE is included in the contact interac- 
tions along with the other short-distance effects. 

We derive the RGEs for the spin-singlet and spin- 
triplet channels from the effective averaged action ansatz 
up to including the 0{p'^) in the expansion of momenta 
and the pion mass. The nontrivial fixed point of phys- 
ical importance is found to be the same as that in the 
pionless NEFT in each channel. The eigenvalues (scaling 
dimensions) and the corresponding eigenoperator of the 
linearized RGEs around the fixed point are also shown 
to be the same. That is, there is one relevant contact 
operator to be resummed. The other operators should 
be treated as perturbations. The L-OPE is also treated 
as a perturbation. A part of the S-OPE contained in the 
relevant operator is resummed to all orders. 

We emphasize that the effects of pions do not alter the 
scaling dimensions and hence the power counting. The 
pions affect the details of the RG flows. We show that in 
the spin-triplet channel the flow is more affected by the 
pions and the region of the initial values that flow to the 
strong coupling phase is larger than in the spin-singlet 
channel. 

We believe that the difference between these channels 
eventually leads to the existence of the bound state (the 
deuteron) in the spin-triplet channel, and the nonexis- 
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tence in the spin-singlet channel. 



B. Comments on the related works 

In the following, we discuss the relation of the present 
paper to the relevant works in the literature. 

As stated in Introduction, our work is very closely re- 
lated to the work by Beane, Kaplan, and Vuorinen [2^ . 
who, working with the PDS renormalization scheme, in- 
troduce a regularization mass scale A to separate the pion 
exchange into its long-distance and short-distance parts, 
and the short-distance part is represented as contact in- 
teractions. The separation of the pion exchange into two 
parts is essential to improve the convergence of the EFT 
expansion in the spin-triplet channel and is similar to 
our Wilsonian RG analysis, with the regularization scale 
A corresponding to our floating cutoff A. 

There are however several points to be addressed: (i) 
Even though they regard A as a low-energy scale, they 
actually consider high-momentum values in their numer- 
ical calculation, ranging from 600 to 1000 MeV. (ii) They 
employ the PDS renormalization scheme simply assum- 
ing the modification does not affect the power counting, 
(iii) The renormalization scale /i and the regularization 
scale A seem to play a similar role in reordering the EFT 
expansion but they are treated independently. As a re- 
sult, A becomes just a new parameter, (iv) They con- 
sider the regularization only for the spin-triplet channel 
because of the singular l/r^ potential, but not for the 
spin-singlet channel. 

From our point of view, these may be seen as follows: 
(i) Our separation scale A is smaller than the physical 
cutoff Ag « 400 MeV, so that it can be consistently re- 
garded as a low-energy scale, (ii) We show in this paper 
that the nontrivial fixed points as well as the scaling di- 
mensions are the same as those in the pionless theory. We 
have shown that the scaling dimensions near the nontriv- 
ial fixed points in the pionless theory correspond to those 
of the PDS ones [s^]. (hi) In our Wilsonian RG analy- 
sis, there is only one scale A, which plays the role of /i, 
renormalization scale, in the RGEs, and simultaneously 
the role of A, the separation scale, (iv) We think that one 
should consider such a regularization for both channels 
from a general EFT point of view. 

The effects of pions in the Wilsonian RG context have 
been studied by Birse [l^ . He shows that the range 
of momenta in which the expansion converges is very 
small in the spin-triplet channel and claims that the pion 
exchanges should be resummed. He then employs the 
distorted- wave RG [1^ and finds that the scaling dimen- 
sions shift by —1/2 (in our convention) in the spin-triplet 
channel at the nontrivial fixed point. Even at the triv- 
ial fixed point, the scaling dimensions are found to shift 
from the canonical values. He also notes that the scale A^ 
(which is the same as Ann in Eq. (|4.10p ) is numerically 
small (^ 2m^) so that the resummation is justified. 

We note that he docs not introduce the separation scale 



A for the pion exchange consistently with the contact 
interactions, but introduces an additional regularization 
scale (i?) and keeps it fixed when studying the RG flows. 
Thus his distorted- wave RGE does not take into account 
the contributions from the OPE to the contact interac- 
tions (diagrams (b) to (d) in Fig. [1]) at all. As a result, 
his RGE does not have a smooth transition to that of 
the pionless theory, where all the pion exchange effects 
are represented as contact interactions. That is, the pion 
exchanges never decouple. In contrast, our RGEs have 
a smooth transition to those of the pionless theory, as 
shown in Appendix \^ and the pions decouple as they 
should. In addition, because of the iterative property of 
the RGEs, the definition of the dimensionless coupling, 
7, is uniquely determined, and it leads to the pertur- 
bative treatment of L-OPE, as explained at the end of 
Sec. HvH 

Fleming et al. [l3| also report the nonconvergence of 
the EFT expansion in the P- waves (the ^Pq channel). Al- 
though we have not explicitly calculated the RGEs in this 
channel, we think that a similar machinery also works 
there. The difference between the S-waves and the P- 
waves is that the physical NN system is near the non- 
trivial fixed point in the S-waves, while in the P-waves 
it seems to be near the trivial fixed point. Although the 
scaling dimensions obtained in the PDS scheme with the 
pole at Z? = 3 subtracted are shown to be the same as 
those at the nontrivial fixed point in the S-waves, they 
correspond neither to those at the nontrivial fixed point 
nor to those at the trivial fixed point in the P-waves [13 . 
Thus the simple PDS with the pole at = 3 subtracted 
should not be used. 



C. Prospects of future research 

Finally we would like to make a comment on a possible 
implementation of the findings of the present paper into a 
more tractable way of calculating the physical amplitudes 
to higher orders. The Wilsonian RG method with the 
momentum cutoff is theoretically transparent but practi- 
cally too complicated to do higher order calculations. A 
simple but powerful scheme that is also consistent with 
our results is desired. Such a scheme would employ the 
dimensional regularization. Since dimensional regular- 
ization does not have a natural separation scale in itself, 
one should introduce it by hand. Thus it would be very 
similar to the BKV prescription. It seems necessary, how- 
ever, to make a connection between the renormalization 
scale fi and the separation scale A so that we have a con- 
sistent RGEs with those obtained in the present paper. 
Work in this direction is now in progress. 

Appendix A: RGEs for the case of A < m,r 

In Sec. lIVI we have derived the RGEs for the case m^r < 
A by expanding the contributions in powers of m^/A. 
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But the pionful NEFT is valid also for the case A < m-^. 
In this Appendix we present the RGEs for the case A < 
TOtt . The diagrams which contribute to the RGEs are the 
same. The difference is that the contributions are now 
expanded in powers of A/m^. 

In the spin-singlet channel, we have 



dx 

dy 
dt 



dz 
'dt 



(first line of Eq. (|TBa|) ) 

— 2(a; + y + z + u)'y — 7^, 
(first line of Eq. (|46b|) ') 

3 

-{2x + 3y + z + 2u)j - -7^, 



(first line of Eq. (j4.6cp ) 

+ (a- + y-z + u)7 + -7 , 



du 



-u — 2{x + y + z)u — 2u7, 

where we have introduced new notations, 
A^ , ml 



7 = —7, u 
mt 



A2 



Now the RGE for 7 is given by 
d^ 



dt 



= -37. 



(Ala) 
(Alb) 

(Ale) 
(Aid) 

(A2) 
(A3) 



Similarly, in the spin-triplet channel, we have 

dx' 
~dt 



dy^ 
dt 



(first line of Eq. (|4.12ap ) 
-2(x' + ^t' + y' + z'-4^(7')7-97^ (A4a) 
(first line of Eq. (|4.12bl) ) 

-{2x'+ 2u'+ 3y'+ z'- Aw')^ - y72,(A4b) 



dz 

— (first fine of Eq. (|4.12cp ) 



+ (.t' + u' + z/ - 2' - 4w')7 4- -7^, (A4c) 



dw' 
dt 

du' 
'dt 



(first line of Eq. (|4.12dp ) 

+ {2x' + 2u + 2y' + 2z' - w')7 + 27^ (A4d) 

-u' - 2{x' + y' + z')u' - 2u'7. (A4e) 



(c) 

In the case of A < TOtt, the operator corresponding D2 
cannot be distinguished with the operator corresponding 

(c) 

to Cq , because these operators are of the same form to 
this order and the pion mass is not a small parameter to 
be expanded. Thus they appear only though the combi- 



nations Cq'^' -|- m'^D'^' , or, ;\; = a; -|- u and x' = x' -\-u' . In 



terms of these variables, the RGEs can be rewritten as 

X) 



dt 



dy 
dt 



dz 



(first hue of Eq. (|4.6ap with x 

-2(x + y + ^)7-7^, 

(first line of Eq. (|4.6b() with x 



(A5a) 



X) 



^(2x + 3y + z)7- ^72, 



(first line of Eq. (j4.6cp with x 
1 



X) 



^)7+ ^7^ 



+{x + y 

for the spin-singlet channel, and 
dt 



(A5b) 



(A5c) 



dji 
dt 



(first fine of Eq. (j4.12ap with x' - 
-2(x' + ^/' + 2'-4w')7-97^ 
(first fine of Eq. (j4.12bp with x' - 

-(2x' + 3y' + z'-4w')7- y7', 



X') 



(A6a) 



X') 



(A6b) 



dz' 
Hi 



— ^ (first fine of Eq. (j4.12cp with x' ^ x') 
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4«/)7+ 27', 



(A6c) 



dw' 
~dt 



Hx+y'- 

(first fine of Eq. (j4.12d|) with ^ x') 
-h(2x' + 2y' + 2z'-^^;')7 + 27^ (A6d) 



for the spin-triplet channel. Here we have now included 
the other terms that we neglected in the ansatz to the 
order 0{p^), such as terms proportional to m^. In a 
similar way, we may consider that x ^-nd x' contain the 
terms proportional to m^, but also terms of all order in 
the expansion in to^ . These x E^nd x' should be compared 
to the couplings x and x' in the pionless NEFT. 

The new coupling 7 is a natural measure of the 
strength of the pion exchange for A < m^, as 7 is for 
TOtt < A. Note that the RGE for the 7, Eq. (|A3p . shows 
that the L-OPE is more irrelevant and the effects of pions 
on the contact interactions thus become negligible very 
rapidly in this region. This implies that one may put 
7 = as a good approximation. In this way, the RGEs 
of the pionful NEFT is smoothly connected to those of 
the pionless NEFT. 



Appendix B: The case of QED 

In this Appendix, we briefly discuss the case of QED in 
the nonrelativistic region (NRQED). As a concrete exam- 
ple, we have a hydrogen atom (or, electron-proton scat- 
tering ) in mind, and we are interested in the low-energy 
region where even the electron behaves as a nonrelativis- 
tic particle. 
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The system consists of a nonrelativistic proton and 
an electron interacting by exchanging photons. The La- 
grangian is similar to that of NEFT with propagating 
pions, Eq. p.Sp . Note that we also include contact inter- 
actions of protons and electrons. By a similar analysis, 
one easily finds that the IR enhancement takes place for 
the proton-electron reducible diagrams, giving the lead- 
ing order contributions to the RGEs. There, (the time- 
time component of) the photon propagator is replaced 
with the instantaneous Coulomb potential. (It is in- 
dependent of the choice of the gauge.) Effectively, the 
RGEs are generated by the averaged action consisting of 
the contact interactions and the instantaneous Coulomb 
interaction, such as 

- Y Jdtjd^x d'y(e^(t, x)<72pit, y))\e^{t, x)a2p{t, y)) 
xC{\x-y\), (Bl) 

where C{r) is the Coulomb potential, 



C{r) 



1 

47rr ' 



(B2) 



similar to that of NEFT, Eq. ([3:22)) . The Coulomb in- 
teraction here should be considered as the long-distance 
part of the Coulomb interaction (L-Coulomb), while the 
short-distance part (S-Coulomb) is included in the con- 
tact interactions. 

The difference between the pionful NEFT and the 
NRQED comes from the following facts: (i) the photon is 
massless while the pion is massive, (ii) the Coulomb in- 
teraction is not a derivative coupling, while the NNtt is, 
and (iii) the electromagnetic coupling is dimensionless 
while the NNtt coupling (5^/2/)^ is dimensionful. 

For the NRQED RGEs, we introduce the dimensionless 
coupling constant jqed-, 



where M is the reduced mass of the proton and the elec- 
tron. This particular A-dependence is required by the 
condition that the RGEs for the contact interactions do 
not contain the explicit A-dependence, as in Sec. IIV Al 
for the pionful NEFT. The RGE for -/qed in the leading 
order is given by 



dt 



~1QED- 



(B4) 



IQED = 



M 



(B3) 



(Note that the usual beta function contribution vanishes, 
de^ /dt = 0, because we are in a region where A is smaller 
than the electron mass.) We see that the L-Coulomb 
interactions become more important in the infrared, in 
strong contrast to the L-OPE. Photons do not decouple 
as they should not. There is no "photonless" QED. 

The contact interactions lead to the similar RGEs to 
those of the NEFT. There is a nontrivial fixed point sim- 
ilar to the one found in the NEFT. Since there seems 
to be no fine-tuning for NRQED, it is natural to con- 
sider that the physical system is close to the trivial fixed 
point. Therefore, the S- Coulomb is irrelevant and should 
be treated as a perturbation. 

In the pionful NEFT case, the L-OPE becomes rapidly 
irrelevant for A < m^r, as is shown in AppcndixlX] On the 
other hand, NRQED does not have such a region where 
the L-Coulomb becomes irrelevant because the photon is 
massless. 

Therefore, we conclude that the L-Coulomb interac- 
tion is nonperturbative, and the contact interactions in- 
cluding the S- Coulomb interaction are pcrturbative. It is 
interesting to compare this with the case of the pionful 
NEFT, where the L-OPE is perturbative while (a part 
of) the S-OPE should be resummed to all orders. 
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